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Abstract 

The semileptonic decays, Be — >Xc{hc) + I + Vi, and the two-body nonlep- 
tonic decays, B^ — ^Xc{hc) + h, (here Xc and he denote (ccpPj]) and (cc["'^Pi]) 
respectively, and h indicates a meson) were computed. Ah of the form 
factors appearing in the relevant weak-current matrix elements with Be as 
its initial state and a P-wave charmonium state as its final state for the 
decays were precisely formulated in terms of two independent overlapping- 
integrations of the wave-functions of Be and the P-wave charmonium and 
with proper kinematics factors being 'accompanied'. We found that the 
decays are quite sizable, so they may be accessible in Run-II at Tevatron 
and in the foreseen future at LHC, particularly, when BTeV and LHCB, 
the special detectors for B-physics, are borne in mind. In addition, we also 
pointed out that the decays Be —> he + ■ ■ ■ may potentially be used as a 
fresh window to look for the he charmonium state, and the cascade decays. 
Be XcpA,2] +1 -\- T^i {Be — > XcP-Pi,2] + h) with one of the radiative decays 
XcP-Pi,2] J I'll) being followed accordingly, may affect the observations of 
Be meson through the decays Be J ftp + l + {Be — > J/jp + h) substantially. 

PACS Numbers: 13.20.He, 13.25.Hw, 14.40.Nd, 14.40.Lb, 12.39.Jh 

I. INTRODUCTION 



The meson Be, being a unique meson, contains two different heavy flavors. It decays 
by one of the two heavy flavors through weak interaction and it happens that the two have 
a comparable possibility each other in magnitude, or by the two heavy flavor annihilation, 
hence, its decay-channels which have a sizable branching ratio, are manifested much richer 
than those of the mesons B^, B^, Bg, D^, D^, Dg etc. Therefore one may study the two 
heavy flavors b, c simultaneously with the meson Be alone, as long as its different weak 
decay channels can be distinguished from each other well. Of all the mesons, in studying 
two heavy flavor b, c simultaneously. Be is unique. 

The meson Be is just discovered very recently. The first positive observation was suc- 
cessful in CDF at Tevatron, Fermilab through the semi-leptonic decays Be — + ^ + ^'/i 
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and the mass tjib, = 6.40 ± 0.39 ± 0.13 GeV, the hfetime tb, = OA6t^o\l ± 0.03 ps etc were 
obtained |T|]. 

Before the observation of CDF, i?c- meson production spectroscopy and various 
decays H^TO| had been widely computed. Now the further experimental studies of the meson 
are planned at Tevatron (in Run II) and at LHC etc. Particularly, in addition to CDF, 
DO, ATLAS and CMS, the detectors BTeV and LHCB are specially designed for B-physics, 
numerous events (more than 10*^ ~ 10^° per year) at these two colliders are expected to be 
recorded so a lot of interesting decay channels of Be will be well-studied experimentally, 
and certain rare processes will become accessible. Therefore, further extensive theoretical 
studies of this meson are freshly motivated. 

The semileptonic decays. Be Xc{hc) + 1 + 1^1, and the two-body nonleptonic decays, 
Be — > Xc{hc) + h, i.e. the decays of the meson B^. to a P-wave chamonium state are certainly 
interesting, but still missing in literature, thus we devote this paper to report our the latest 
computations on them, although the semileptonic decays were reported shortly ]10[- Why 
the decays interest people, let us outline the reasons below. 

First of all, people would like to know how sizable the decays will be, especially, to 
know if accessible in Run-II of Tevatron and/or in LHC. Especially the cascade decays of 
Be ^ Xc-\ — ■ and Xc ^ J/i' + 1 looks quite like as a signal for the observation of the meson 
Be through the decays Be J/ip + ■ ■ ■, because the photon may be missed in detectors. 
In addition, two of the P-wave charmonia have a branching ratio about a few tenth for the 
radiative decays XcfPi] ^ J/'ip + l {Br = 27.3%) and XcP^2] ^ J/ip + ^f {Br = 13.5%), 
so indeed the cascade decays may potentially contribute a substantial background for the 
observation Be meson through Be ^ J/ip + ■ ■ ■■ Therefore, even only from the point of view 
to estimate the background for the observation on Be meson, to see how great the concerned 
decays is very interesting. 

If one would like to see CP violations in Be decays, for example, to see CP violation in 
the decays Be —>■ h + hi + h2 {h, hi, h2 denote various possible mesons), as emphasized in 
Ref . |TT[ , one knows that the interference of the direct decays with a cascade one through 
a resonance, e.g., XcP-Po], i-e.. Be — >■ XcP-Po] + h and XcP-Po] hi + h2, may enhance the 
visible CP violation effects substantially. Thus to see the advantage of this method for the 
purpose quantitatively, the knowledge on the decay Be — > XcP-Po] + h is necessary. 

QCD-inspired potential model works very well for nonrelativistic double-heavy systems. 
The systems {cb) and {cb) in forming bound states, except the reduce mass, are similar to 
the well-studied systems {bb) and (cc), so it is believed that with potential model the static 
properties of the systems {cb) and {cb) can be predicted very well as those of bottomium {bb) 
and charmonium (cc). In general, to apply the wave functions to computing the relevant 
decay matrix elements is attracting, since the potential model will have further tests. Thus 
with the wave functions of Be (the ground state of the system of {cb)) and Xc{hc) (the P-wave 
states of (cc)) obtained by potential model, we have applied the wave functions to compute 
the decays Be — >■ Xc{he) + ■ ■ ■■ 

Since the mass of Be {ttlbJ is much greater than those of the P-wave charmonia {m,^^ 
and rrihj, so the momentum recoil in the concerned decays can be a great (even relativistic). 
If one tries to apply the Schrodinger wave functions of nonrelativistic binding systems to 
computing the decay processes with such a great (even relativistic) recoil momentum, one 
cannot carry out the computation of the decay matrix element successfully just as done in 
atomic and nuclear decays by taking a suitable 'reference frame' and then a simple 'boosting', 
since the recoil in an atomic or nuclear process is much smaller than that in the present 
concerned decays. The great momentum recoil obviously means the velocity between the 
two CMS of Be meson and the charmonium state is huge, and the potential wave functions of 
the parent and the daughter states, given just in each CMS respectively, cannot be applied 
directly just by choosing a suitable reference frame and simple boosting the wave functions 
to the same frame. Thus when applying the wave functions to calculation of the decays (e.g. 
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the semileptonic decays and most two-body nonleptonic decays here) with such a great (even 
relativistic) momentum recoil, special handling is needed. To deal with the momentum recoil 
properly, an approach for the decays from a nonrelativistic S-wave state to another 5'-wave 
one, the so-called generahzed instantaneous approximation, was proposed in Ref. [Q]. Since 
it is straightforward to extend from a nonrelativistic S-wave state to another S'-wave one, to 
the present case, that the decays are from a nonrelativistic S'-wave state to a P-wave one for 
the approach, hence here we do so. The key points of the approach may be outlined as the 
three steps: firstly, to 'extend' the potential model, which is based on Schrodinger equation, 
to the one on Bethe-Salpeter (B.S.) equation^ even for the non-relativistic binding systems; 
then, according to Mandelstam method ||12|| to formulate the (weak) current matrix element 
(an elementary factor for the relevant decays) sandwiched by the B.S. wave functions of 
the two bound-state, so that the current matrix element is written in a fully relativistic 
formulation; finally, by making the so-called 'generalized instantaneous approximation' on 
the fully relativistic matrix element i.e. to integrate out the 'time' component of the relative 
momentum in the Mandelstam formulation by a contour integration, and as the final result, 
the current matrix element turns out back to be formulated in terms of proper operators 
sandwiched by the Schrodinger wave functions of the 'original' potential model. Since the 
weak current matrix (by means of the Mandelstam method) was formulated relativistically , 
so we can be sure that the final formulation takes the recoil effects into account properly and 
no new free parameter is added at all. Besides the great recoil effects are treated properly, 
one additional advantage of the approach is that it has a more solid ground on quantum 
field theory than that on the 'original' potential models, because the B.S. wave functions 
and the Mandelstam formulation have a more solid 'ground' on quantum field theories and 
they are used as a starting point to make the generalized instantaneous approximation. 

On the other hand, B.S. equation is four- dimensional in space-time to describe a bound 
state problem, and there are a few problems still, such as, how to determine the QCD- 
inspired four-dimensional interaction kernel of the equation properly, and what is the physics 
meaning of the excitation in its relative-time 'freedom' of the two components etc. In 
addition, the B.S. equation is harder than a Schrodinger one to solve, even when the four- 
dimentional kernel is fixed. Whereas with the generalized instantaneous approximation, 
the current matrix elements are reduced into certain proper operators sandwiched by the 
potential model Schrodinger wave functions finally, therefore, the approach, in the meantime 
to circle the difficulty about treating the great momentum recoil effects properly, has also 
kept some of the advantages of potential model, such as to avoid the difficulty to solve the 
B.S. equations etc. 

Finally we should note here that in our calculating the two-body nonleptonic decays of 
Be to the P— wave Xc and he states, the so-called factorization assumption and the effective 
Lagrangian for four fermions in which the 'short-distance' QCD corrections have been taken 
into account with OPE (operator product expansion) and RGM (the renormalization group 
method), as done by most authors, are adopted. 

The paper is organized as follows: To follow the Introduction in Section-II, the exclusive 
semileptonic differential decay rates, the matrix elements and form factors etc are described. 
In Section III, the adopted approach, the so-called generalized instantaneous approximation, 
to compute the form factors is illustrated precisely. In Section IV, the two-body non-leptonic 
decays of are formulated with necessary description. Finally in Section V, numerical 
results and discussions are presented. The dependence of the current matrix elements on 
the form factors, and the dependence of the form factors on and ^2, the integrations of 
the wave function overlapping, are put in Appendix. 



^For the binding systems, Be and Xc{hc), to do the extension is just by means of the original 
instantaneous approximation pr opo sed by Salpeter, that can be found in many text book on quan- 
tum field theory e.g. the book iQ to 'build' the relation between the Schrodinger equations and 
the relevant B.S. ones. 
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II. THE EXCLUSIVE SEMILEPTONIC DECAYS AND RELEVANT CURRENT 

MATRIX ELEMENTS 

The T— matrix element for the semileptonic decays Be X^c + + I'f. 



Gf _ 



(1) 



where X^s denotes Xc and h^, Vij is the Cabibbo-Kobayashi-Maskawa(CKM) matrix element 
and is the charged current responsible for the decays, p, p' are the momenta of initial 
state Br and final state Xrr- Thus we have: 



(2) 



where hfj^^, is the hadronic tensor and l'^'^ the leptonic tensor. The later Z^,^ is easy to compute 
whereas in general the former /i^^, can be written as: 

V = -(^9,iu + f3++{p + p')^l{p + p')u + + p)^,{p - p')v + P-+{p - p'),,{p + p'),y+ 

/3.4p - p')^ip - p'), + z7W<x(P + P'YiP - P'T , (3) 
and by a straightforward calculation, the differential decay-rate is obtained accordingly: 



dxdy 



2 

mf , 



4M'2 - m2 



4M2 



M2 



r/2 



9M' — 

+4(/3.-+/3-.)^(2-4a; + ,-^^ 



1/(1 



/2 



M2 



, m2 M'2 



(4) 



where x = Ei/M and y = {p — p')^/M^, M is the mass of i?c meson, M' is the mass of final 
state Xcc- The coefficient functions a, 7 can be formulated in terms of form factors. 
Note here that we have kept the mass of the lepton m/ precisely that is different from those 
by N. Isgur et al |P and by B. Grinstein et al [|I4|, so the formula here can be applied not 
only to the cases of e and /x semileptonic decays but also to those of r-semileptonic decays. 
1. If Xcc is hc{[^Pi]) state: the vector current matrix element 



< Xccip', e)\V^\Bc{p) >= re; + s+{e* ■p){p + p')^ + s„(e* ■ p){p - p')^, 
and the axial vector current matrix element 

< Xccip', e)\A^\Bcip) >= tve^,p^e*''ip + p'fi^p - p'Y , 



(5) 



(6) 



where p and p' are the momenta of Be and he respectively, e is the polarization vector of he- 
2. If Xcc is Xc(P-Po]) state: the vector matrix element vanishes, and the axial vector 
current 
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< X,,{p')\A^\B,{p) >=u+{p + p')^ + u4p-p')^. (7) 

3. IfX^gis Xc(PA]) state: 

< X,,{p', e)\V,mp) k; + c+(e* ■p){p + p'), + c_(e* ■ p){p - p'), , (8) 

and 

< e)\A^\B,{p) >= tqe^^p^e^^ip + p'Y{p - p'f . (9) 

4. IfX,g is Xc(pP2]) state: 

< X,,{p',t)\V^\B,{p) >= ih+-t^,p„e*^"p^{p + py{p-pr , (10) 

and 

< x,,{p', e)\A^\B,{p) >^ ke;y + b+{e;yp'^){p + p'), + b4e;yp'^){p - p')^. (ii) 

The form factors r, s+, s_, f , m^, m_, /, c+, c_, /c, and are functions of the mo- 

mentum transfer t = {p — p'Y can be calculated precisely. In Ref. |^ we proposed an 
approach, the generalized instantaneous approximation, to compute those form factors for 
the decays of B^. to an S'-wave charmonium state J /ip or r/c- Now we are computing the form 
factors r, s+, s_, ■ ■ • appearing in the decays of Be to a P-wave charmonium state, in fact, 
the approach may be used directly, thus the approach is adopted in the present calculations 
here. 



III. THE SO-CALLED GENERALIZED INSTANTANEOUS APPROACH TO THE 

WEAK CURRENT MATRIX ELEMENTS 

To calculate these form factors, the approach developed in Ref. |§] is adopted. Let 
us outline the approach here for convenience. According to the Mandelstam formalism 
when the considered weak (electromagnetic) current matrix element involves only one 
hadron in the initial state and one in final state respectively, then it may be written down in 
terms of Bethe-Salpeter (B.S.) wave functions which describe the hadrons as bound states 
exactly: 



(27r)4 



Xp'(g')r'^Xp(9)(^2 + m2)l , (12) 



where Xpio)-, Xp'W) ^ire the B.S. wave functions of the initial and final states with the 
corresponding momenta p, p' . Throughout the paper we use pi, p2 denote the momenta of 
the quarks in the initial meson S^, and p'^, p'^ denote the momenta of the quarks in the final 
meson or h^. For convenience let us introduce further definition of the relative momentum 



q \px q'): 



mi 

Pi = aip + q, ai 



P2 = a2P - q, a2 



nil + 

1712 



mi + m2 

Pi,P2,mi and m2 are the momenta and masses for the quark and antiquark respectively. 
Note that the matrix element of the current Eq.(12) now is fully relativistic, thus it can be 
used as the start 'point' to take into account the recoil effects in the decays no matter how 
great the recoil moment is in the considered decay. To prepare in applying the generalized 
instantaneous approach for the matrix element, we need to 'convert' the potential model 
onto the B.S. equation 'ground'. 
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A. The Potential Model and B.S. Equation 

In general, the B.S. equation for the corresponding wave function Xpil)'- 

r d^k 

{A - mi)xp{q){h + 1712) J J^:^^(P^ ^' (l)Xp{k), (13) 

where V{p, k, q) is the kernel between the quarks in the bound state, may describe the 
relevant quark- antiquark bound state well. Accordingly the B.S. wave function Xp{q) should 
satisfy the normalization condition: 



d'^qd'^q' 



■Tr 



Xp{q)^^ [S^\pi)S2\p2)5\q - q') + Vip, q, q')] Xp{q') 



= 2ipo , (14) 



where Siipi) and 82(^)2) are the propagators of the relevant particles with masses mi and 
1712 respectively. 

As pointed out in introduction, the B.S. equation in four dimension should be reduced 
to a one in three dimension i.e. the time-like component momentum should be integrated 
out (the instantaneous approximation) with a contour integration as proposed by Salpeter, 
especially when the kernel has the property as follows 

Vip,k,q)^V{\k -q \) , 

to do it is very easy. When one make a contour integration of the 'time' component of 
the relative momentum on the whole B.S. equation, then the B.S. equation is deduced 
straightforwardly into a three-dimensional equation which just is a Schodinger equation 
in momentum representation. Since the start point of the common potential model is a 
Schodinger equation, thus we may convert the potential model onto a ground based on the 
B.S. equation in the way with instantaneous approach. 

To treat the possible great recoil effects in the decays, furthermore we need to convert the 
instantaneous approximation to a covariant way too, i.e. to divide the relative momentum q 
into two parts, q\\ and q±, a parallel (time-like) part and an orthogonal one to p, respectively: 

q' = <ii + <± , 

where g^y = ip-q/Mp)p^ and q^j^ = q^ — q^^y Correspondingly, we have two Lorentz invariant 
variables: 

p-_q 

In the rest frame of the initial meson, i.e., P= 0, they turn back to the usual component 

go and I 9 I, respectively. 

Now the volume element of the relative momentum h can be written in an invariant 
form: 

d^k — dkpkpj^dkpTdsdcl), (15) 

where is the azimuthal angle, s — {kpqp — k ■ q)/ (kpTqpr) ■ Now the interaction kernel can 
be rewritten as: 

V{\k-q\)^V{kp^,s,qp^). (16) 

Defining: 

'Ppiqp±) = ^ j ^Xp{qp\\,qp±), 



(ipT^^q^p-q^^^f^^- 



kprpdkpxds 



^P^' ~ J (27r)2 
The B.S. equation now can be rewritten as: 

Xp{qph(lp±) = Si{pi)r]{qp±)S2{p2) 
and the propagators can be decomposed as 



Siipi 



+ 



J{i)qp + UiM - ujip + ie J{i)qp + aiM + ujip - ie ' 



with 



ip 



m' 



+ <lpT,^ip{<lp±) 



2lj. 



ip 



M 



Uip ± J{i){mi + 



lp±) 



(17) 



(18) 



(19) 



(20) 



where i = 1,2 and J{i) = (—1)*+^. Here A^(gpj_) satisfies the relations 



0. (21) 



Due to these equations, may be referred as the p— projection operators, while in the rest 
frame of corresponding meson, they turn to the energy projection operator. 
We define V3p^(q'p_L) as 



V^J^(gp±) = ^%{.qpi.)j^ 



'^pi.(lpi-)T7^2p 



M 



(22) 



where the upper index C denotes the charge conjugation. In our notation, {qpi) = 
A2p{qp±). Integrating over on both sides of Eq . ([T8|) , we obtain: 

(M - uip - uJ2p)<^^^{qp±) = Atp{qpi.)Vp{(lpi.)^2p{(lp±y^ 

(M + uip + uJ2p)fp~{(lp±) = ^ipiqp±)Vp{(lpL)^tpi(lp±y^ 



<^p (gp±) = Vp^{qp±) = 0. 

The normalization condition of Eq. (p^ ) now becomes 



(23) 



tr 



2Pn . 



From these equations, one may see that in the weak binding case to compare with the 
factor [M — ijj\p — uj2p)-, the factor (M + a;ip + a;2p) is large, so the negative energy components 
of the wave functions <y9 are small. In the present case, for the heavy quarkonium and 
meson, it is just the case, so we ignore the negative energy components of the wave functions 
safely at the lowest order approximation. 

Neglecting the negative energy components of the wave functions, the B.S. equation 
contains the positive component 

^p^iqp±) = A^p(gp±)— ^p(gp±)— A2jf (gpx) 
7 



only, and the normalization condition becomes: 



tr 



_++ i) 

^ M 



2Po 



Now let us consider the wave function appearing in the above equations. We know 
that the total angular momentum of a meson is composed from orbital one and spin, fur- 
thermore there are two ways i.e. S-L coupling or j-j coupling to compose the total angular 
momentum. Here to consider P-wave states of charmonium, we adopt the way of S-L cou- 
pling, i.e. we let the spins of the two quarks couple into a total spin, which can be either 
singlet or triplet, then the total spin couple to the relative orbital angular momentum, and 
finally we obtain the total angular momentum. In this way, the reduced B.S. wave function 
(/7p can be written approximately as: 



2V2M 



for ^5*0 state, and 



(24) 



(25) 



for ^5"! state, where e"^ is the polarization of this state. For the P-wave (cc) wave functions: 

^-pM) = ^^y^-f5^nlMM, (26) 

for ^Pi, i.e. he state, and 

P -^ A/f . 

(27) 



->■ P + M 

"p'pM) = TTT^ ^(WniM.(?) < IS,,LM,\JJ, >, 



2^/2M 

for ^Pj{J=0, 1, 2) i.e. Xc states, where e is the polarization vector of total spin, 
< ISz, LMz\JJz > is Clebsch- Gordon coefficients which couple L, S to the total angular 
momentum J. tpnoo and ipniM^ are the full B.S. wave functions. 



B. The Radius B.S. Equation in Momentum Space 

To solve the B.S. equation, the key problem is about its radial component. If we ignore 
the negative energy contributions, the reduced B.S. equation Eq.(18) in the rest frame of 
the meson center mass system can be written as: 

= M-u,-u, ■ 

In the frame, the energy projection operator: 

A+ = -^(a;i7o+ 7 ■ ? +mi), 
= 7^('^27o- 7 • ?-m2). 
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where the kernel V acts on ip{(l) as: 



(29) 



i.e. to correspond to the potential model more precisely, the interaction kernel can be 
formally divided into the corresponding non-perturbative QCD 'linear' one, Vs (in scalar 
nature) and the corresponding gluon exchange one, (in vector nature). 

When substituting Eqs. (^ , |26 ), the wave functions in the meson center mass system, to 
the reduced B.S. equation Eq.(^), the equation for a spin singlet state 5 = becomes: 



05=o(5) 



AuJiUJ2{M — iUi — iU2) 



mim2 



4V,{q,k)-W,{q,k) 



^s=o{k)dk\, (30) 



where the (j)s=o{q) is (l)nooC'So) or <^niA/^(^-Pi)- Since square of the relative momentum q is 
small to compare with quark mass squared in the 'double heavy' meson, as a lowest order 
approximation, we have ignored such higher terms and use uJi = mi, uj2 = ^2 in numerator. 

Now let us factorize out the radial component of the wave function and its relevant B.S. 
equation in momentum space from the angular ones: 

i^nLMA^) = <PnL{\ Q |)lLAf, (^, V') , 

where n is the principal quantum number, L is the orbital angular momentum and is the 

projection of the third component of L, 0nL(| ^ |) is the radial wave function and Ylm^O, 0) 
is the spherical harmonic function. For the spin singlet states, multiplying Ylf^^q) to two 
sides of the reduced B.S. equation and sum over by using the formula, 

4:71 



where 6 is the angular between the unit vector q and k, the radial reduced B.S. equation for 




^Sn state is obtained: 



0no(| 9 I) 



AuJiUJ2{M — iOi — UJ2) I 
Whereas for ^Pi state: 

1 



|mim2 J 



WJq,k)-mq,k] 



M k\)dk\. (31) 



Ml 1 I) 



AuJiUJ2{M — LUi — UJ2) 



mim2 



AVAq,k)-Ws{q,k] 



k \)cos9d k 



(32) 



where (/)no(| 1 I) and (j)ni{\ |) are the radial parts of the wave functions. 
Similarly, for the spin triplet states 5* = 1 we have: 



1 



Im 



Im 



AuJiUJ2{M — iUi — UJ2) 



X < mim2 



AV,{q,k)-Ws{q,k] 



<Ps=i{k)d k 



(33) 



where the 



^5=1 



q) is 



^Si) or 



^Pj). Then the equation for ^Si is: 
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and for ^Pj: 

M\Q\) = 



4:UJiUJ2{M — uji — UJ2) 



AuJiUJ2{M — uJi — UJ2) 



|mim2 J 



mim2 



4V,{q , k) - Ws{q , k] 



0„o(| k\)dk\; (34) 



</>ni(| k \)cos6d k 



(35) 



The normalization of (pnL now is read: 

q^dqj 



17111712 
UJ1UJ2 



2M. 



Under the present further approximation, the three triplet P-wave states ^Pj and the 
singlet ^Pi as well, are degenerated. The reason is that we have ignored the 'splitting' 
interactions at all. 



C. The Generalized Instantaneous Approximation 



After neglecting the negative energy component and the 'treatment' above, the weak 
current matrix elements become as follows: 



d^q 



Tr 



0J\ 



le 



qp + aiM — 001 + ie 



xv{(lp±)- 



A^(gp±) 



-qp + a2M - UJ2 + ie_ 



(36) 



The generalized instantaneous approximation, being an extension for the original one on 
the B.S. equations suggested by Salpeter, with the Cauchy's theorem performs a contour 
integration about the time-like component qp in complex plan on the whole current matrix 
elements precisely. As the final result, the matrix elements turn out to become a three 
dimensional integration about the space-like components q±. 

If we choose the contour along the lower half plane, after completing the contour inte- 
gration, the current matrix elements become as follows: 



(Pqi 

{2^y 



-Tr 



vWp'±)^iWp' 
M' -uj[ 



-r 



^, ^^{qp±)v{qpi-)^t{qpl^ 

M — uji — UJ2 



This matrix elements can also be written in the frame where the momentum q'^ is the 
integral argument by means of a suitable Jacobi transformation, i.e. 



Qp'Tdgp'T^S 

(2vr)2 



tr 



^'^ilp'±)^t^^t^i(lp±)j^ 



(37) 



The above formula with the argument q'j_ as the integral argument is more convenient, 
especially, in the cases when we calculate the matrix elements involving a P-wave state in 
the final state. 

After performing the calculations on the matrix elements precisely, the dependence of 
the matrix elements on the overlapping integrations of the initial and the final state wave 
functions becomes transparent. So is all the form factors too. 

Since there is the so-called spin symmetry for heavy mesons, all of the form factors for 
their decays may attributed to one 'universal' function i.e. the Isgur-Wise function [jl5 
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Therefore for the double heavy meson Be to a 5*- wave charmonium, at the hmiting >> 
rric >> Aqcd i-G- turning to the case of the heavy mesons, the form factors are attributed to 
the Isgur-Wise function, and the Isgur-Wise function is related to an overlapping integration 
of the wave functions of Be and the S'-wave charmonium with certain kinematics factors 
precisely [P|,p!B|. Now at the present case of B^ to a P-wave charmonium, not only due to 



the spin-symmetry but also due to the great recoil in the decays, alternatively there are two 
independent and 'universal' functions, essentially, just two overlapping integrations of the 
wave functions of the initial and final bound states, and and all of the form factors 
are described by the two general functions with proper kinematics factors precisely. 

Since in the present case the initial state is of an S-wave and the final state is of a P-wave, 
so the matrix elements must be related to two kinds of terms: one is to the integration which 
does not depend on the relative momentum qp,j_ at all, and the one just on g^,^ linearly. 
Namely all the form factors appearing in the decays depend on two universal functions 
and ^2 only: 



(27r)= 



etm^ = / ^^€iMM'T)i^noo{qpT)q';^, (38) 



where 



{p-p'r 

M'2 



M2 



describes the polarization vector along recoil momentum p', e^{L) is the polarization vector 
of the orbital angular momentum. 

We should note that for the decays from an S-wave state to a P-state, the function 
generated in the present approach is special. Since has more direct roots to the momentum 
recoil, so it cannot be obtained by boosting the final state wave function as done in the cases 
with a small recoil. The reason is that approaches to zero when the momentum recoil 
vanishes. Whereas, the function as in the cases with a small recoil, can essentially involve 
recoil effects just by 'boosting' the final state wave function. 

Substituting the B.S. wave functions Eq.(p4D and Eqs.( ]26| - p7D into the equation of current 



matrix elements and using Eq.(^), the precise formula for the form factors i.e. the precise 
dependence of the form factors on ^1 and ^2, can be obtained and we put them in the appen- 
dices and the curves of ^1 and ^2 obtained by numerical calculations are shown in a figure. 
With the functions .^i, .^2 and the form factors, the decay rates of the semileptonic decays 
and the spectrum of the charged lepton for the decays can be obtained by straightforward 
numerical calculations. 

Note that in our calculation on the form factors, we have used the relations: 



Y.ei{S)et'{L)<lX;lX'\00 >= .^{g^ 

\,X' » 



A,A' 



, 2M' 

E <{S)e':{L) < lA; iy|2A" ^^J), (39) 

A,A' 
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where < 15*2; lLz\JJz > as previous are C.-G. coefficients. The polarization vector 
J), J = 1 and the tensor e^^{J), J = 2 have the projection properties: 

A 

E = liKc^Kf^ + KpPlc) - \pluPLp- (40) 



IV. THE TWO-BODY NON-LEPTONIC DECAYS 

In this section we outhne how the two-body non-leptonic decays B^. — * Xc{hc) + h (here 
h denotes a meson) are calculated. We adopt the factorization assumption on the decay 
amplitudes which is widely adopted in estimation of the non-leptonic decays for various 
mesons. With the assumption, the weak current matrix elements appear in the calculations 
precisely and they are related to the form factors just obtained in the previous section. For 
the non-leptonic decay modes — > Xc{hc) + h (caused by the decay b ^ c), the following 
effective Lagrangian Lg// (QCD corrections are involved) is responsible: 

L^fj = ^ { V4cMQf + C2ifi)Qf] + h.c.} 

+ penguin operators . (41) 

Gp is the Fermi constant, Vij are CKM matrix elements and Cj(/i) are scale-dependent 
Wilson coefficients. The four-quark operators Ql'^ and (CKM favoured only) are: 

Qf = [V:, {du)v^A + C {su)y.A + idc)v-A + (scW-a] (c6)v-^ 
Q2 = Kd {cu)v-A {db)v-A + {cu)v-A {Sb)v-A + Krf {cc)v-A (db) 

+ {cc)y^A m, (42) 

where (gig2)y-A denotes gi7^(l -75)'?2- 

Because at this moment we restrict ourselves to consider the decays in which the coef- 
ficients of 'penguin' operators in the effective Lagraingen are small in comparison with the 
two main ones Ci and C2, so the contribution from penguin terms is neglected in the calcula- 



tions, although in the Ref. fl^ it is pointed out that in total decay width the penguin may 
have interference with the main ones and can course an increase about %3 ~ 4. Moreover, 
at this stage we also restrict ourselves only to consider the decay modes where the weak 
annihilation contribution is small due to precise reasons e.g. the helicity suppression etc0, 
namely we neither take into account the contribution from the weak annihilation here. 

Precisely by means of the factorization assumption, the decay amplitudes for the 
non-leptonic decays can be formulated into the three factors: the so-called leptonic de- 
cay constants, which are defined by the matrix elements: < 0\Afj,\M{p) >= ifMPfj. (or 
< 0|V^|V^(p, e) >= fyMye^)] the weak current matrix elements < Xc|V^(^/x)|-Sc >, which 
are those as the semileptonic decays; and the relevant coefficients in the combinations: 
ai = Ci + KC2 and 02 = C2 + kci, here k = l/N^ and A^^^ is number of color. The coefficients in 



^We will consider the contribution from penguin and weak annihilation carefully elsewhere. 



12 



the combination ai, 02 is due to the weak currents being 'Fierz-reordered'. In the numerical 
calculation later on, we will choose oi = Ci and 02 = C2, i.e., we take k = in the spirit of 
the large Nc limit, and QCD correction coefficients ci and C2 are computed at the energy 
scale of nib- 

Therefore with the relations between the currents and form factors obtained as in 
the semileptonic decays, finally the factorized amplitudes for the nonleptonic decays can 
be formulated in terms of the form factors and the decay constants by definitions: < 
0\A^\M{p) >= ifMPfi and < 0| V"^|V^(p, e) >= fyMyefj,. Thus the decay widths for the 
two-body nonleptonic decays can be computed straightforward. 



V. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we present the numerical results. 

In the numerical calculations, based on potential models the parameters are chosen as 
follows: 

A = 0.24 GeV2, a = 0.06 GeV, Aqcd = 0.18 GeV, a = e = 2.7183, Vo = -0.93 GeV, 



Vbc = 0.04 PI, mi = 1.846GeV, m2 = 5.243GeV. 



With this set of parameters, we obtain the masses: 

Mb, = 6.33 GeV, M' = 3.50 GeV, 

and corresponding radial wave-functions of Be meson and P-wave charmonium Xc he nu- 
merically. Here in the present evaluations, we only carry out the lowest order ones without 
considering the splitting caused hj L — S and S — S couplings, in which all the bound states 
^Pj{J = 0, 1, 2) and ^Pi are degenerated. 

To see the behaviors of the universal function ^i{tm — t) and ^2{tm — t) i.e. the two over- 
lapping integrations of the wave functions of initial and final states, we plot them explicitly 
in Fig.l, where = (M - M')^ t = {P - P'f. 

The lepton energy spectra for the decays i?c — >■ Xc + + i^, for which the mass of 
charged lepton can be ignored, are shown in Fig. 2, and those for the decays ^ Xc + i' + 
for which the mass of charged lepton r cannot be ignored, are shown in Fig.3, where \pe\ is 
the momentum of lepton. The difference between Fig. 2 and Fig.3 is due to the sizable mass 
of r-lepton. For the semileptonic decays, we put the corresponding widths in Table I. 

As for the non-leptonic two-body decays Be Xcihc) + h, we only evaluate some typical 
channels, whose widths are relatively larger, and put results in Table II. In the numerical 
calculations, we have chosen ai = ci and 02 = C2, i.e., k = 0, and ci and C2 are computed at 
the energy scale of m^. The values of the decay constants: = 0.131 GeV, fp+ = 0.208 
GeV, fa, = 0.229 GeV, fK+ = 0.159 GeV, fK*+ = 0.214 GeV, fo, = 0.213 GeV, /dj = 0.242 
GeV, //)+ = 0.209 GeV, = 0.237 GeV are adopted by fitting decays of B and D mesons. 

if comparing the results in Table 1 with the decays of Be to 5'-wave charmonium states 
J/ip and rje e.g. T{Be ^ J/ip + 1 + u) ^ 25 ■ 10~^^GeV [H,^, one can realize the semileptonic 
decays of Be to the P-wave charmonium states in magnitude are about tenth of the decay 
Be ^ J /ip + I + vi. As for the two-body nonleptonic decays, due to the difference in 
momentum recoil and the fact that the recoil momentum is fixed in a given specific decay, the 
P-wave decay Be — >■ Xc{he) + h can be greater than twentieth of the one. Be — > J/ip{rie) + h, 
to an S-wave state. 

The ffist observation of Be by CDF group is through the semileptonic decay Be 
J/ijj + I + ui, hence, we can conclude that most of the decays concerned here are accessible 
in Run-II of Tevatron and in LHC, especially, when the particular detector for B physics 
BTeV and LHCB at the two colliders are concerned. It is because that Tevatron and LHC 
will have more than 20 time events of Be meson than Run-I and have much better detectors. 

Since the decays Be Xc['^-Pi,2] + I + i^i have such a quite sizable branching ratio, so 
the cascade decays i.e. the decays with an according one of the radiative decays XcP-Pi,2] 
J/ip+'J followed may affect the observation through the semileptonic decays Be J /ip+l+Pi 
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as done by CDF group substantially, especially, when the efficiency of detecting a photons 
for the detector is not great enough. 

We also would hke to point out here that with sizable branching ratio, the decays Be — > 
he + I + and/or Be he + h potentially can open a fresh 'window' to observe the 
charmonium state hd^Pi], especially, to note that the charmonium state hd^Pi] has not 
been well-estabhshed experimentally yet. 
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APPENDIX A: 

In this appendix, we present the form factors and formulas for a, f3^^ and 7 which are 
required in the calculations on the exclusive semileptonic decays of Be to Xce, which denotes 
one of ^Pi, ^Pq, ^Pi and ^P2 states as indicated precisely in each case below. 

For convenience, we introduce the parameters below: 

_ , p-p' 
L1J20 = LO2' 



'MM' 



uio = \luJ2Q -ml + ml , 



1. Be Meson to Charmonium hd^Pi] 

The matrix elements for the vector and axial currents: 

< X{p',€)\V^\Be{p) >= re; + s+{e* ■ p){p + p')^ + s_(e* ■p){p-p')^ 

< X{p',e)\A^\Be{p) >= ive^.pae*''{p + p'np-pT- 

Where 

(m[ - m2){mi + ujiq -1712- 0^20)6 

r = 

8m'ia;ioa;2o 

_ jm'i + m2){mi + ujiq + m2 + a;2o)6(p • p') 
d>M'Mm'^ujwOJ2Q 

_ m2[M{m2 + UJ2Q - mi- ujiq) - M'{mi + uj^q + m2 + a;2o)]6 
m2[M(m2 + UJ20 - nil - ^10) - M'{mi + cuio + m2 + uj2o)]i\ 



d>M' MuJiQUJ2onep 

m2[M{m2UJ2o + '^20 ~ 'm,iUJ2o - ujIq) - M'{miUJ2o - ujIq + m2a;2o + '^2o)]^2 



(Al) 



8M'M^u;fo^^2o 
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+ 



{m[ + 7712) (mi + cjio + 1712+ a;2o)6 i^' - M)^i 



16M'Mm'iu;iou;2o SM'M^uio ' 

m2[-M{m2 + 0020 - rui - ojw) - M'{mi + uio + 1712 + a;2o)]6 

m2[—M{m2 + 0020 — mi— uio) — M'{mi + coio + m2 + Ci^2o)]Ci 

8M' Mu!ioU!2o'nep 



m2[-M{m20J2o + ojIq - miijj2o - u;fo) ~ M'{mxUJ2Q - ool^ + m2a;2o + ^2o)]6 

SM'M'^uIqU20 

(m'l + 1712) (mi + uJiQ + m2 + 0^20)6 , - 



+ 



8M'M2u;io ' 



_ (m'l + m2) (mi + m2 + a;io + c<;2o)6 
~ IQM' Mm'iUJwUJ2Q 

The dependence of a, and 7 on the above form factors: 



/3. 



M2 



(i-y)-i 



I 71J-2 p 2 



"'2 /i#/2a„,2 1 



M2 



M'2 



:i-z/)-3 



1 

+ 4 



M'2 



P 



4M'2 



+ 



M2 



M'2 



:i-y)-3 



I 71^2 2 

rs- + M j^s_ 



7 = . 
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2. Be Meson to Charmonium XcP^o] 

The matrix element for the vector current vanishes in the present decay. 
The matrix element for the axial current: 

< X{p')\A^\B,{p) >= u+{p + p')^ + u4p- p')^. 

Where 

[M{ujiQ + mi + U20 + m2)Ci + M'{uj2o + m2 — o^io - rni)]m2^i 



8y/3M'uJioUJ2onep 
M'(a;io + mi + 0020 + m2)Ci + M{u;2o + m2 — coio — mi)^i 



SVSM'uionep 

3^2 [^'(^'1 + m2){uJio + mi+ 0020 + ^2) + M{m[ - m2)(cu2o + m2 - ooio - mi)] 

16^/3M'Mm'lCc;loa;2o 

_^^2'm2[M' {-ujio - mi + U20 + m2) + M{u2o + m2 + uJio + mi)] 

8^/3M'Mu}ioU}^o 

^ ^2[M'm2{m2U}2o + a;|o - mia;2o - cj^q) + Mu2oim2U}2o + oojo + mia;2o - oufo)] 

8^/3M'MoufoU}2o 

^^2[-M'{m2u;2o + ^20 + "^i'^20 - i^^o) + M{-m2UJ2Q - 0^20 + '^1^20 + <^fo)] 

^ UM'-M){m2 + u;2o) 
8V^M'Mm2(jJio ' 

_ [-M{oJio + mi+ 0020 + "^■2)6 + M'{uj2o + m2 - ooio - mi)]m2^i 

8\/3M'uJioUJ2onep 
M'{ujio + mi+ 0020 + ^2)6 - M{uj2o + m2 - (jJw - mi)^i 

SVSM'uoionep 

3^2[M'{m[ + m2)(a;io + mi + UJ20 + "^2) - M{m[ — m2){u;2o + m2 — coio — mi)] 

16V3M'Mm[u;iou;2o 

^2m2[M' (-ojio - mi + 0020 + ^2) - M(u;2o + m2 + o^io + mi)] 

^^2[M'm2{m2UJ2o + - ''^i'^2o - t^?o) ~ Mcj2o("^2'^2o + '^io + "^1*^20 - ^?o)] 

8\/3M'Mo;?oCc'2o 

^2[-^'("T'2'^20 + (^lo + "^1*^20 - t^?o) ~ M(-m2a;20 - ^^20 + "^1*^20 + 'AoA 



%\PiM'Mujl 
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^ 6(M^ + M)(m2+a;2o) .^^3) 
8V3M'Mm2C<;io 



The dependence of a, and 7 on the above form factors: 
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q; = 0, (A13) 

/3-+ = U-U+, 13— = ut , (A14) 
7 = 0. (A15) 

3. Be Meson to Charmonium XcpA] 

The matrix elements for the vector and axial current currents: 

< X{p',e)\V^\B,{p) >= le; + c+{e*-p){p + p')^ + c.{e*-p){p-p')^, 

< X{p', e) \A^mp) >= iqe,,p,e*\p + p'Y{p - p'Y . 

Where 

^ ^ (mi + uw + m2 + c^2o)6[(P • P'f - M'^M''^]m2 
4:V2MM'^nepuioU2o 
(mi + ujio + 1712+ a;2o)6(p • P') ("^1 + <^io - m2 - a;2o)6(p • P') 



2\^MM'u}ioU}2o 'i\/2Mm\m2UJiQ 



Ay/2m'iUJwUJ2Q Ay/2M'^M'^uj 



10 



^ i2m2[{p ■ p'f - MH'I'^] 

4\/2M2M'2cUiOCU20 



(cuio + mi + CJ20 + m2) (miCc;2o + m2C^2o + ^20 ~ ^10) 

' 2 ' 

(A16) 



_ (mi + uw + m2 + uj2o)^i{M''^ - P ■ p')m2 ^ (M'^ - p ■ p')^2 
^+ Tn^TTTTTo ^ 



^ (mi + uio + m2 + a;2o)C2 _^ (^i + a;io - m2 - 0^20)6 



10 



a;ioc<;2o 8-v^Mm'im2Ci;io 



^2m2{M''-p-p') 
8V2M'^M'^uJwUJ2o 



{uio + mi + U20 + ^2) (mia;2o + m2a;2o + i^jo - (^ip) 

2 

<^20 "^lO 



^ ^ (mi+a;io + m2+a;2o)Ci(^'^+P-p0^2 ^ (M^^+p-p0^2 
{nil + <^io + ^2 + uj2o)C2 {mi + ujio - m2 - 0^20)6 



4:\^MM' 00100020 8V2Mm[m200io 



^ C2m2{M'^ + p ■ p') 
8V2M^M'^oowo02o 



{ooio + mi + 0020 + 1712) (mia;2o + m2a;2o + 0020 - ooIq) 

2 

^20 ^10 
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(tuio + mi + UJ20 + ^2)6 'm2{-uJio -1711+ ^20 + "^2)6 
S\/2M'uJionep 8\/2M'uJioUJ2onep 

{LJ20 + miUJ20 + m2UJ20 - (^fo)^^ , ("^2 + c^2o)6 



8V2M'M(jj^ 



+ 



10 



8V2M'Mm2(jJw 



6^i2 



8V2MM'a;ioa;2o 



(-a;io - mi + W20 + m2) , (-miu;2o + m2Cc;2o + 1^20 - ^10) 



^20 



+ 



10 



The dependence of a, and 7 on the above form factors: 



(A17) 



M2 



M2 



M'2 

-►2 



(i-y)-3 



1 

+ 4 



(AI8) 



(A19) 



(A20) 



4M'2 



M2 



M'2 



:i-y)-3 



7 = 2/g . 



(A22) 



4. 5c Meson to Charmonium XcP^2] 

The matrix element of the vector and axial currents: 

< X{p', e)\V^\B,{p) >= ih+_e^,p,e*''y^{p + p'y{p - pl)\ 



< X{p', e)\A^\B,{p) >= keiy + 6+(e>V)(p + + &-(e>''p")(P - 
Where 

^ _ (mi + uJiQ + m2 + oJ2o)ii ^ m2(- mi - low + m2 + <x;2o)6 



4a;ionep 



A.ujiQUJ2Qne-P 
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6"^2 



4:Mujiouj2o 



-ujio - mi + U20 + m2) 

^20 



+ 



AMm2U0iQ 
-miUJ2Q + m2UJ2Q + ujIq - ujIq) 



10 



_ m2{mi + Uio + 1712 + UJ2q)^1 
0+ — TTTTTT H 



6 



8M' MuJiQUJ2onep 



6^2 
8M2M'a;ioa;2o 



(wio + mi + c^2o + ^2) , (mia;2o + m2a;2o + 0J20 - ^^w) 



^20 



+ 



UJ 



10 



6_ = 



6 



6"^2 



8M2M'cuioti^2o 



m2(mi + cjio + m2 + a;2o)Ci 
8M' Mu;iQUJ2onep 

{ujio + mi + UJ20 + m2) , (mia;2o + m2a;2o + a;|o - uj^q) 



+ 



^20 



^10 



m2(mi + c^io + m2 + cj2o)6 
8M' MuJioUJ2onep 



6 



6"^2 



8M2M'a;ioa;2o 



(a;io + mi + LJ20 + 1712) {miU}20 + m2UJ2o + <^io - <^fo) 



^20 



+ 



The dependence of a, and 7 on the above form factors: 



M2 p' 
2M'2 



(A23) 



Ap' 



26^ p' 



\ 



+ 



^2 



(A24) 



A/f2 „' i-2 



1 - 



^2\ ^4 

M2 4M2p' I 26+6_M>' 



V 



M'2 M'4 



M'4 



+ 



M2p 
6M'2 



^2 

6M^2 



(A25) 



+- 



(A26) 



M2j)' ,2 



2M'2 



2(M^ + M'^) - M^y 



2hi_ p' 
~3 



+ 
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3M'2 



M2 
M'2 



(l-y) + 3 



(A27) 



7 



M^p' kh 



(A28) 
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TABLES 



TABLE I. The semileptonic decay widths (in the unit 10~^^ GeV) 







T{Bc-^^ Police) 




r{Bc^''P2iue) 


e(M) 


2.509 


1.686 


2.206 


2.732 


r 


0.356 


0.249 


0.346 


0.422 



TABLE II. Two-body non-leptonic B+ decay widths in unit lO'^^ GeV 



Channel 
-ft)" 


r 

a? 0.569 
„2 fi q-|7 

a? 0.0815 
af 0.277 


r(ai = 1.132) 
0.729 

0.104 
0.355 


Channel 
^Pip 


r 

a? 1.40 

til u.ouu 
af 0.331 
af 0.579 


r(ai = 1.132) 

1.79 

1 .uo 
0.425 
0.742 




a? 1.71 
a| 1.03 
af 0.671 
al 1.05 


2.19 
1.33 
0.859 
1.34 


3PnK+ 
3PiK+ 
3P2K+ 


a? 4.26 X 10-3 
a\ 2.35 X 10-3 
af 0.583 X 10-3 
af 1.99 X 10-3 


5.46 X 10-3 
3.02 X 10-3 
0.747 X 10-3 
2.56 X 10-3 


^PiK* 
^PqK* 
^PiK* 
^P^K* 


a? 7.63 X 10-3 
af 4.43 X 10-3 
af 2.05 X 10-3 
af 3.48 X 10-3 


9.78 X 10-3 
5.68 X 10-3 
2.63 X 10-3 
4.47 X 10-3 


^PlDs 
^P2Ds 


a? 2.32 
af 1.18 
a? 0.149 
af 0.507 


2.98 
1.51 
0.191 
0.650 


^PiD* 
^PoD* 
^PiD* 
^PiD* 


a? 1.99 
af 1.48 
af 2.21 
af 2.68 


2.56 
1.89 
2.83 
3.44 


^PlD+ 
3PoD+ 
3PiP)+ 
3P2L>+ 


a? 0.0868 
af 0.0443 
a? 0.00610 
af 0.0209 


0.111 
0.0568 
0.00782 
0.0267 


^PiD*+ 
^PoD*+ 
^PiD*+ 


a? 0.0788 
af 0.0567 
af 0.0767 
af 0.0972 


0.101 
0.0726 
0.0983 
0.124 
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FIGURES 




1 2 3 4 5 6 7 



FIG. 1. The universal functions and ^2 

vs. — t. They are the overlapping-integrations of 
the the wave functions for Xc(^c) and Be with the definition as in Eq.(38). The sohd hne is of ^1, 
the dashed one is of ^2- 




0.5 1 1.5 



IPelGeV 



FIG. 2. The energy spectrum of the charged lepton for the decays ^ Xc + + where 
the solid line is the result of /ic[^-fi] state, dotted-blank-dashed line is of XcP-Fb], dashed line is of 
XcP-Pi], dotted-dashed line is XcP-P2]- 
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FIG. 3. The energy spectrum of the charged lepton for the decays Be Xc + t + Vr, where 
the soUd line is the result of hd^Pi] state, dotted-blank-dashed hne is of XcP-Fb], dashed hne is of 
XcP-Pi], dotted-dashed line is XcP-fb]- 
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